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Abstract 


Soil and water conservation is a priority in the watersheds of arid and 


semi-arid regions for the proper planning and integrated management of 
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water resources. The objective of this work is prioritizing 91 watersheds 
in 14 regions of Peru with different geographical, hydrological and 
geological conditions, through the integration of Geographic Information 
Systems (GIS), Principal Component Analysis (PCA) and the Weighted 
Sum Approach (WSA). In addition, homogeneous regions were identified 
with hierarchical Cluster Analysis in R and Ward's method. The result 
showed the existence of 19.49 % of total area in high and very high 
priority category in two regions. The results of the Cluster Analysis 
showed that 35 % of the watersheds have homogeneous zones within the 
geometric and shape factor, while 65 % within the drainage and relief 
factor. In general, GIS, PCA and SWA methods are an efficient tool, which 
allows decision-making authorities for better planning and conservation 


of natural resources in hydrographic watersheds. 


Keywords: Erosion, morphometric, drainage, cluster analysis, factor 
analysis. 


Resumen 


La conservación del suelo y el agua es una prioridad en las cuencas 
hidrográficas de las regiones áridas y semiáridas para la adecuada 
planificación y gestión integrada de los recursos hídricos. El objetivo de 
este trabajo fue priorizar 91 cuencas hidrográficas en 14 regiones del Perú 
con diferentes condiciones geográficas, hidrológicas y geológicas, 
mediante la integración de sistemas de información geográfica (SIG), 
análisis de componentes principales (PCA) y el análisis de suma 
ponderada (WSA). Se identificaron regiones homogéneas con el análisis 
de conglomerados jerárquico en el método de R y Ward. El resultado 
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mostró la existencia de un 19.49 % del área total en categoría de alta y 
muy alta prioridad en dos regiones. Los resultados del análisis de 
conglomerados mostraron que el 35 % de las cuencas hidrográficas tiene 
zonas homogéneas dentro del factor geométrico y de forma, mientras que 
el 65 % dentro del factor de drenaje y relieve. En general, los métodos 
SIG, PCA y SWA son herramientas eficientes, que permiten a las 
autoridades encargadas la toma de decisiones para una mejor 
planificación y conservación de los recursos naturales en las cuencas 


hidrográficas. 


Palabras clave: erosión, morfometría, drenaje, análisis de 


conglomerados, análisis factorial. 
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Introduction 


Peru, globally, is one among the countries most affected by variability and 
global climate change (Vega, Lavado, 8l Felipe, 2018) impacting on the 
spatial distribution of water resources in Peru's Andes, a fundamental 
source of water supply for the population, agriculture and energy 
production (Wongchuig, Mello, 8 Chou, 2018). In this context, there is a 
requirement for sustainable watershed management (Bhattacharya, 
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Chatterjee, € Das, 2019) ranging from an efficient morphometric 
characterization (Ghosh 8 Gope, 2021) with new scientific tools to enable 
adaptation to variability and global climate change with attention on 


prediction until the solution of problems (Jacobs 4 Brian, 2020). 


The Principal Component Analysis (PCA) method allows within the 
various morphometric variables to explain the most significant part of the 
variance of the data (Helness, Damman, Sivertsen, € Ugarelli, 2019) to 
identify cause - effect relationships; zones with homogeneous climate or 
similar hydrological characteristics. This analytical approach groups 
similar data while  ¡identifying relationships between variables 
(Gorgoglione, Gioia, € lacobellis, 2019) precisely to identify the most 
dominant ones in hydrological processes (Balbín et a/., 2020). Subyani, 
Qari and Matsah (2012) performed the PCA and Cluster Analysis (CA) but 
only to ten watersheds of different sizes in the western region of Saudi 
Arabia with 18 morphometric variables. Yunus, Oguchi and Hayakawa 
(2014) presented the geomorphological quantification of mountainous 


terrain, but limited to 36 watersheds in the Western Arabian Peninsula. 


Research such as that of Gajbhiye and Sharma (2017) is based on 
the prioritization of watersheds based on the PCA and the Composite 
Factor (CF), this includes the calculation of the simple arithmetic average 
of the final priority rankings, but they have the disadvantage that they 
assume the same degree of importance to all morphometric parameters 
(Aher, Adinarayana, € Gorantiwar, 2014; Bharath et a/., 2021). The 
aforementioned aspect is contradictory because each watershed has its 
own characteristics of shape, relief and hydrographic network. Recent 
research by Malik et a/. (2019), and Setiawan and Nandini (2021) have 
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shown strategic optimism between PCA and WSA techniques for 
watershed prioritization. The WSA method is known because it offers 
consistency in land surface analysis that contributes to watershed 
prioritization. According to Singh and Singh (2018) emphasize that WSA 
is a traditional watershed prioritization method compared to the rest and 
provides dynamic and efficient results. New integration methods are 
essential that allow better management of hydrographic watersheds in 
Peru, due to its climatic diversity and fragile ecosystems. Therefore, the 
objective of the research is to prioritize watersheds for soil and water 
conservation based on GIS, PCA and WSA techniques in different regions 


of Peru. 


Materials and methods 


Fourteen regions with different hydrographic, hydrological, geological and 
ecological conditions were selected as the study area for the 91 
watersheds. These watersheds are located in two of the country's largest 
watersheds, the Pacific and the Atlantic. The digitized watersheds were 
processed in the GIS of the ArcMap 10.5 software, and based on the 
formulas in Table 1, the morphometric parameters of the watersheds can 
be calculated. The characterization of the watersheds will be carried out 
in qualitative and quantitative terms that are grouped into surface, shape, 


drainage and relief parameters. 
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Table 1. Formulas of the morphometric parameters of the watersheds. 


Parameter Equation Units References 
watersheds amplitude W = AL¿* km Horton (1932) 
Slope of the watersheds Se = e» NA % Zavoianu (1985) 

Gravelius compactness 
K.= 0.282P4”"* - Zavoianu (1985) 
coefficient 
Form factor F¿ = AL - Horton (1932) 
Shape index Su == LA % Horton (1945) 
Elongation ratio Ro =11284A "17" - Schumm (1956) 
Miller, Ritter and Kochel 
Circularity ratio R¿ = 4mAP"? - 
(1990) 
Relief ratio Ry = HL - Schumm (1956) 
Mean stream slope S= AHI” % Zavoianu (1985) 
Drainage density Di =L¿¿7* km/km2 Horton (1945) 
Mean length of overland 
Le=0:5D3* km Horton (1945) 
flow 
Constant channel 
C =1Dz* km/km2 Schumm (1956) 
maintenance 
Roughness coefficient E, = DaS; - Horton (1945) 


Open Access bajo la licencia CC BY-NC-SA 4.0 
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Principal component analysis 


It is a difficult and complicated process to classify potential areas with 
more benefits from water resources in 91 watersheds. The aim of principal 
component analysis (PCA) is to explain the variance-covariance structure 
in multiple data sets using a few linear combinations of the original 
variables, according to Kottegoda and Rosso (2008). The PCA technique 
is designed to transform p correlated X variables, which are known, into 
an equal number of uncorrelated Z indices. PCA is a technique for 
summarizing a complex set of data, distinguishing quantitative dependent 
and independent variables, as well as identifying closely related variables. 
Matrix notation in its general form is represented using Equation (1). In 


which Z and X are nxp matrices and A is a pxp matrix of coefficients: 


Where Z;is an nx1 (column) vector and a; is a px1 (column) vector 


of coefficients., j of 1, 2, ... p. 


Considering that each column is an observed variable in rigor of n 
answers to the same question, one can study the variance and covariance 
between different variables, represented by a square C matrix of order 
pxp symmetric (Yang et a/., 2020). The variance-covariance matrix X ¡is 
represented using equations (2) and (3). In practice is estimated by the 


sample covariance matrix C: 
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C11  C12 Cip 

o A C2zp 
> : 3 (2) 

Cp1 Cpz Cop 

1 

Ci = 721 Xi Xi (3) 


Where C becomes the sample correlation matrix; X” is the transpose 
of the standardized X variables, and c;; is the value of each element of the 


covariance. 


In the PCA, the variance of the Principal Component (PC) is 
maximized, which is mainly related to the eigenvalues of the sample 
covariance or correlation matrix (Malik et a/., 2019; Kottegoda € Rosso 
2008). The correlation r between the X¡; of the PC is obtained by the ratio 
of the covariance of the X matrix and the variance of Z; of the PCs using 
Equation (4): 


Cov(x;¡;zj) 


r(x;52;) > JVar(z)) (4) 


Where r is the correlation coefficient; Cov (xX;, z;) is the covariance 
of X variables in the ¡th row and jth row of the PCs; and, Var (z;) is the 
variance of jth PCs. Value of r greater than 0.90 indicates strong 


correlation, r between 0.75 and 0.90 indicates good correlation, and r 
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between 0.60 and 0.75 indicates moderate correlation (Sharma, 
Gajbhiye, € Tignath, 2015). 


Factor analysis 


Factor analysis (FA) is a procedure similar to principal components 
analysis to identify physically significant variables and reduce 
morphometric parameters. The analysis via this technique produces easily 
interpretable results. FA techniques are well described in the literature 
(Finch, 2020). The correlation between the variables is related to the 
factor loadings, considering correspondence between the factor score 
matrix and the Zs in the PCs, the orthogonal in the rotation has an impact 


on the PCs and the corresponding factors are equivalent. 


Therefore, a comparison with the unrotated components is possible, 
while the rotation should be such that meaningful physical interpretations 
are possible from the resulting components obtained (Jolliffe, 2002). 


Thus, the rotated components are given by using Equation (5): 
Z¿=XS  RPDAOS (5) 


Where Z; is the rotated components; D-%3 is a diagonal matrix in 
which the nonzero elements are the square roots of the eigenvalues of 
the correlation matrix S expressed by Equation (2), and RP the orthogonal 


rotation. 
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Cluster analysis 


The hierarchical Cluster Analysis (CA) technique was applied to classify 
the homogeneous groups called clusters. The method applied for data 
agglomeration was in R mode and the Ward's criterion (Ward, 1963), one 
of the most widely used in hydrology and meteorology (Cupak, Walega, 
e: Bogustaw, 2017), with the objective of achieving similarity between 
watersheds based on morphometric parameters with a focus on water 


resource and soil. 


Weighted sum approach 


WSA was applied to the most significant morphometric variables that are 
identified with PCA, WSMP or morphometric parameter weight was used 
in cross-correlation analysis (Siddiqui, Said, € Shakeel, 2020) whereas, 
Composite Factor (CF) was evaluated for final priority ranking and related 
category, which are expressed according to Singh and Singh (2018), Malik 
et al. (2019), and Setiawan and Nandini (2021) as: 


_ Sumof correlation coef ficient 
WS COAH 


(6) 


Total of correlation 


CF = PRoempWsmp (7) 
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Where CF is the composite factor; PRsmp is the preliminary ranking 
of the most significant morphometric parameter of the PCA, and Wesmp is 


the Weight of the significant morphometric parameter. 


Based on the CF value, the watershed was assigned a priority 
ranking. The lowest CF value belongs to priority ranking one, the second 
lowest value to priority ranking belongs to priority ranking two and so on 
for all watersheds. The study used the SPSS 26.0 statistical package and 
STATISTICA 12 as the computational tool for the PCA. 


Results and discussion 


The results of the statistical analysis in Table 2 of the morphometric shape 
parameters such as compactness coefficient, shape factor, elongation 
ratio and circularity factor indicate the existence of elongated and circular 
watersheds for their mean values between 0.31 to 1.42, with round oval 
watersheds prevailing according to the criteria by Horton (1932) and 
Schumm (1956). In the geological aspect, the circularity factor ranged 
from 0.24 to 0.73, the high results are indicative of rock structures that 


are impermeable and non-homogeneous and that control drainage. 
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Table 2. Statistical analysis of morphometric parameters. 


Parameters Min Max Mean Sd Cv 
Zmax 124.45 5 500.00 3 698.89 1724.57 46.62 
Zmin 4.89 4 013.00 1 516.20 1 183.44 78.05 

H 45.56 5 429.55 2 182.67 1.395.439 63.93 
A 45.90 1 496.85 459.69 372.18 80.96 
P 33.02 258.12 101.92 52.75 51.75 
Lo 20.40 268.38 85.02 12.60 35.98 
L, 49.53 1467.02 449.55 364.45 81.07 
Sc 6.90 50.99 23.34 13.70 58.69 
wW 2.22 21.89 11.30 3.07 49.26 
Kc A, 17 2.06 1.42 0.18 12.49 
Ff 0.10 0.46 0.31 0.09 28.32 
Sw 2.17 9.90 3.68 1.74 47.47 
Re 0.36 0.77 0.62 0.01 15.81 
Re 0.24 0.73 0.51 0.11 21.09 
Ra 0.82 189.18 66.99 44.38 66.25 
S;- 2.36 37.31 11:03 8.90 70.33 
Da 0.98 3.05 1.88 0.39 20.81 
Ls 0.16 0.51 0.28 0.07 23.78 
C 0.33 1.02 0.56 0.13 2346 
C, 1.49 143.78 23.91 25.89 108.28 


Open Access bajo la licencia CC BY-NC-SA 4.0 
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The rugged relief is evidenced by average slopes of 23.34 %, 
classified as mountainous watersheds prone to dominant environmental 
processes such as deterioration of soil fertility, erosion and sediment 
transport in watercourses according to Kumar, Das and Das (2020). 
According to Amare, Kassie and Sulla (2020) significant variations in relief 
influence the distribution of rainfall and other climatic factors. The 
drainage density in the hydrographic watersheds ranges from 0.98 
km/km2 to 3.05 km/km? for an average of 1.88 km/km?. According to the 
categorization of Horton (1945), the watersheds have a medium drainage 
trend and consequently a medium permeability and infiltration of water; 
as well as, moderate tendencies to avenues that induce floods. Results 
that are in correspondence with those obtained by Mahala (2019), but in 
the Kosi River watersheds of North India for the mountainous flat tropical 


environment. 


In general, these are watersheds with a greater probability of using 
water resources spatially and temporally; therefore, there is a high need 
for hydraulic works to mitigate the problems of flooding and frequent 
landslides in the watersheds of Peru's territory. Table 3 and Table 4 show 
the results of the Preliminary Ranking (PRsmp) based on morphometric 


parameters by region. 
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Table 3. Preliminary priority ranking based on watersheds. 


LI 
U r 
> > o == 0) 
y w 
1 A All ale PF iw|53|S 
2 y | 31£8|9 ce ao c|s|k| 9|0 3 
3 N o =, 3 ny 5 Ss|3Y|3 0 pa = | 
o 2 [an |3[|9|0 < ce [> 2/3 | S=|0|9|5y 
e 3 > y = 010 la] ee (2) O pl = 
o (7) y a Y >) 
y 
Zmax 2 13 7 9 11 8 14 6 4 10 | 12 1 3 5 
Zmin 3 9 11 13 14 [10| 4 7 2 8 12 1 5 6 
H 2 13 1 7 5 9 14 6 10/|12 | 11 4 3 8 


Lr 2 8 13 7 12 | 5 4 (14 | 3 1 9 4 6 6 


W 1 5 12 8 14 | 3|11 13| 4 2 9 6 7 10 


Sc 11 9 8 2 7 [12 14/13 3 1 5 4 6 10 


Kc 4 4 12| 9 3 8 5 |10| 2 7 (13| 1 6 11 


Sw 14 11 4 7 11 |13| 3 6 8 [12/10 2 9 5 


Re 11 10 2 7 12|6 9 3 ¡13| 5 1 14 8 4 


Ls 6 3 4 12 1 2 |14 11|9 8 7 5 13 | 10 


53 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 14(6), 40-75. DOI: 10.24850/j-tyca-14-06-02 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


o) 1) Check for updates 
OPEN ACCESS 
Tecnología y 


CienciaszAgua 


Table 4. Preliminary priority ranking based on regions. 


Regions CP PRsmp 
Amazonas 5.53 2 
Ancash 7.79 8 
Apurimac 8.37 10 
Cajamarca 7.89 9 
Cusco 8.74 12 
Huancavelica 6.89 5 
Huanuco 9.21 14 
Junin 8.68 11 
La Libertad 6.37 4 
Lima 5.42 1 
Loreto 8.79 13 
Pasco 5.89 3 
San Martin 6.58 6 
Ucayali 7.78 7 


The PCA was started by obtaining the matrix of the simple linear 
Pearson correlation coefficient (r), there is a strong correlation (r > 0.9) 
between the morphometric parameters of 1) A, P; A, Lc; A, Lr and A, W; 
2) P, Le; P, Lr and P, W; 3) L¿, Lr and Le, W; 4) L, and W; 5) K¿ and Rc; 6) 
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Fr, Sw and Fr, Re; 7) Sw; Re; 8) Da, Lr and Da, C; and 9) Lf and Cc. While, 
there is a good correlation (0.75 < r <0.9) between 1) Hm, Ra and 2) S;, 
Cr. In addition, there are moderately acceptable Pearson correlation 
coefficients (0.60 < r < 0.75) at 1) Zmax, Zmin; Zmax, Ra; 2) P, Kc; P, Rc; 3) 
W, Fr; W, Sw and W, Re. The parameters of S¿ and S, did not show a 
relationship with the rest to the other variables that explain their effect, 
with except for C,. Similar results achieved Yunus et al. (2014) in 
watersheds of the Western Arabian Peninsula, and Mokarram and 
Sathyamoorthy (2015) in hydrographic watersheds of Iran fundamentally 
in the parameters of A, Lc, Re, Fr, Re and Da with strong correlation. While, 
with a certain approximation by those of Malik et a/. (2019) strongly 
correlated to Lr, Fr, Re, Kc, Re and differs with Da in Indian watersheds. 
Based on the significant correlations, no drainage correlations are 
observed with the surface variables, which are in agreement with those 
reached by Sharma et al. (2015). Other authors, such as Gede, Anwar 
and Lasminto (2017), demonstrated that A and L, influence the maximum 
flow and the unit hydrograph at its peak time with a correlation coefficient 


of 0.98, an aspect that could be related in the current research. 


The application of the FA to the twenty morphometric parameters 
in the first iteration, six PCs were defined with eigenvalues greater than 
one and that represent around 88.26 % of the total variance. However, 
the morphometric parameters with weak loads (S. and S,) were extracted 
to improve the FA. Under these conditions, the same factor load was 
reached with 6 PC that explain 94.58 % of the total variance. The Varimax 
rotation showed improvement in the loads, Table 6, which is because the 


matrix of the rotated PC was considered. 


55 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 14(6), 40-75. DOI: 10.24850/j-tyca-14-06-02 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


o) 1) Check for updates 
OPEN (3 Y ACCESS 
Tecnología y 


CienciaszAgua 


In Table 5, the factor load of PC1 constitutes 26.56 % of the total 
variance in the rotated matrix, the second PC2 17.87 % and the third PC3 
explains 16.82 % of the total variance, having these first three 
components a strong charge of 60.95 %; while, from PC4 to PC6, factor 


loadings only explain 33.63 % of the total variance. 


Table 5. Factor loading matrix of reduced morphometric variables 


(Rotated Varimax). 


Principal components 

Parameters A 
PCi1 PC2 PC3 PC4 PC5 PC6 
Zmax 0.15 0.09 -0.10 0.75 -0.03 0.63 
Zmin -0.01 0.04 -0,13 -0.05 -0.07 0.98 
H 0.18 0.08 -0.01 0.97 0.02 -0.05 
A 0.93 0.26 -0.03 0.01 0.21 0.01 
P 0.83 0.32 -0.02 0.02 0.43 -0.01 
Lo 0.94 0.18 -0.03 0.01 0.21 0.01 
L- 0.93 0.27 -0,04 0.01 0.20 -0.01 
wW 0.83 0.52 0.01 0.02 0.17 0.01 

Kc 0.32 0.01 -0.05 | -0.01 | 0.93 | -0.06 | 
Fr 0.28 0.94 0.06 0.03 0.01 0.02 
Sw -0.26 -0.93 -0.05 -0.06 -0.01 -0.05 
Re 0.29 0.95 0.06 0.03 0.01 0.02 
Re -0.28 0.01 0.04 0.07 -0.92 0.03 
Ra -0.28 -0.02 0.03 0.93 -0.09 -0.04 
Da 0.02 -0.08 -0.97 -0.01 0.08 0.07 
Lr 0.01 0.04 0.99 -0.03 -0.01 -0.05 
C -0.01 0.03 0.98 -0.03 -0.01 -0.05 

C, -0.50 0.04 -0.15 | 0.34 | 0.28 | -0.14 | 
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The factor load of PC1 in Table 5 is strongly correlated with the 
factor loadings of A, P, Lc, Lr and W, which could be called geometric 
factors that explain 26.56 % of the total variance of the 18 parameters 
analyzed. The second PC2 has a strong correlation with the variables of 
Ff, Sw and Re, which is why it is called the form factor of the hydrographic 
watersheds that explains 17.87 % of the total variance; while, the PC3 
load factor, its strong correlation is with the Dg, Lf and C called drainage 
factor that explains 16.52% of the total variance. The fourth factor load 
of PC4 has a strong correlation with the variables Of Zmax and Hm, which 
is why it is called the relief factor of the watersheds that explains 13.82 
% of the total variance. The first four components account for 74.77 % of 
the total variance, which explains the most important morphometric 
parameters and the most significant physical factors. Being very useful, 
the PCA to define the most relevant variables. Yunus et a/. (2014) 
achieved results similar to the current ones, in 36 watersheds in which 
three CPs explained 73 % of the total variance that strongly reflected the 
dimensions and surface of the watershed, as well as the drainage texture. 
In the same sense, Sharma et al. (2015) three PCs defined the drainage, 
slope and shape; later with similar results it was by Malik et al. (2019) 


with strong correlation in shape and drainage parameters. 


Table 6 shows that there are variables with a strong correlation in 
each PC. These parameters are considered the most significant for 
applying the WSA model and in turn the prioritization of the watersheds. 
In the final prioritization of the watersheds by region, the composite value 


was calculated based on the preliminary ranking and weight of the six 
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variables Lc, Re, Lr, H, Kc and Zmin with the cross-correlation analysis Table 


6. The composite value was determined using Equation (8). 
CF = 0.26L¿ + 0.23R¿ + 0.10Lf + 0.16H + 0.18K, + 0.08Z min (8) 


Table 6. The cross-correlation between the important parameters. 


Parameter Lo Re Lf H Kc Z min 
Lc 1.00 0.43 -0.01 0.17 0.48 0.02 

Re 0.43 1.00 0.10 0.16 0.11 0.04 

L; -0.01 0.10 1.00 -0.03 -0.07 -0.17 

H 0:17 0.16 -0.03 1.00 0.09 -0.11 

Kc 0.48 0.11 -0.07 0.09 1.00 -0.13 
Zmin 0.02 0.04 -0.17 -0.11 -0.13 1.00 
Sum 2.08 1.84 0.82 1.28 1.47 0.63 
Sum total 8.12 8.12 8.12 8.12 8.12 8.12 
WSA 0.26 0,23 0.10 0.16 0.18 0.08 
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Prioritization of watershed using PCA-WSA 


The final priority ranking of the watersheds was according to the lowest 
CF value as shown in Table 7 and Table 8, where five priority category 
rankings are established. The Amazonas region with very high priority 
category of < 3.10 and area of 11.12 %; the Lima region with high priority 
category of 3.11 < 5.10 but with the minimum affected area of 8.37 %; 
the regions of La Libertad, Huancavelica, Pasco, Ancash, San Martin with 
medium priority category 5.11 < 7.10 and area of 29.71 %; the regions 
of Apurimac, Junin, Ucayali, Huanuco, Loreto with very low priority 
category and the highest percentage of 39.35 % in the least affected. 


Table 7. Final priority ranking of watersheds based on CF value in 


regions. 
Regions CF Final priority ranking 
Amazonas 3.07 1 
Ancash 6.22 6 
Apurimac 9.15 10 
Cajamarca 8.62 9 
Cusco 8.28 
Huancavelica 5.58 
Huanuco 9.83 13 
Junin 9.33 11 
La Libertad 5.57 3 
Lima 5.07 2 
Loreto 10.45 14 
Pasco 5.59 
San Martin 7.09 
Ucayali 9.55 12 
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Table 8. The priority category for watersheds of the regions. 


DS Priority . Area 

No. Priority level Regions 
category (%) 
1 <3.10 Very high priority Amazonas 11.12 
2 3.11 to 5.10 High priority Lima 8.37 


La Libertad, Huancavelica, 
3 5.11 to 7.10 Medium priority 29.71 
Pasco, Ancash, San Martin 


4 7.11 to 9.10 Low priority Cusco, Cajamarca 11.45 


Apurimac, Junin, Ucayali, 
5 > 9.10 Very Low priority 39.35 
Huanuco, Loreto 


The priority ranking map for the watersheds by regions of Peru ¡is 
shown in Figure 1. In the research the highest category level represents 
the highest degree of erosion potential and potential areas for the 
application of soil and water conservation, as well as a need for immediate 
reforestation. In this sense, the Lima and Amazonas regions are in the 
high and very high priority category with a totally vulnerable area of 19.49 
%. According to Puno and Puno, (2019), and Setiawan and Nandini, 
(2021), structural soil and water conservation techniques should be 
applied in the most vulnerable watersheds to minimize the susceptibility 


of the watersheds. 
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Figure 1. Priority ranking map for watersheds of regions in Peru. 


(080 61 


AM 2023, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 14(6), 40-75. DOI: 10.24850/j-tyca-14-06-02 
(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


o) 1) Check for updates 
OPEN ACCESS 


Tecnología y 


CienciaszAgua 


The medium and low priority category is indicative of a moderate 
degree of erosion and partial existence of vegetation cover, in seven 
regions of Peru: La Libertad, Huancavelica, Pasco, Ancash, San Martin, 
Cusco and Cajamarca, representing an area of 41.16 %. Reforestation 
programs and better sustainable use of soil and water resources should 


be implemented in these watersheds. 


In the very low priority category, it is indicated that the watersheds 
have sustainable morphometric characteristics such as the regions of 
Apurimac, Junin, Ucayali, Huanuco and Loreto, in the latter the protection 


of the vegetation cover of the soil should be maintained. 


Homogeneous hydrological characteristics 


The Figure 2 shows the dendrogram obtained from the cluster analysis in 
R mode. All parameters of the watershed are presented with the exclusion 
of light loads (S. and S;,). It is observed that there are two main clusters 
(A and B) and each one is subdivided into two main clusters. The cluster 
A is subdivided into cluster A(I) with two watersheds parameters (Ff and 
Re), and cluster A (II) includes six parameters (A, Ly, Lc, P, W and Kc); 
while, cluster B is subdivided into cluster B (III) that includes two 
parameters (Ls and C) and cluster B(IV) that contains eight parameters 
(Zmax, Zmin, Hm, Rec, Ra, Sw, Da and Cy), the latter, with the lowest 
correlation and highest distance. In this sense, cluster A presents better 


correlation and less distance below 0.2 with respect to cluster B. 
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According to Eltahan, Elhamid and Abdelaziz (2021), although FA is not a 


grouping method, there are similarities between main cluster A (1 and II) 
and CP (1 and 2), as well as B (Il and IT) and CP (3, 4, 5 and 6). 


Ward Linkage: Correlation Coefficient Distance 
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Figure 2. Dendrogram of cluster analysis results in R-mode. 
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Figure 3 shows the result of the dendrogram obtained by the CA for 
the 91 hydrographic watersheds in Peru, numbered in the abscissa axis. 
It is observed that there are two main clusters C and D. Cluster C is 
subdivided into C-1 and C-II with 19 and 13 watersheds respectively. The 
cluster D also subdivided to cluster D-I11 and D-IV in watersheds 11 and 
48 respectively. From the CA of Figure 2 and Figure 3, it can be seen that 
35 % of the hydrographic watersheds have homogeneous areas in 
geometric and shape factor, of which 72 % belong to the Atlantic slope. 
However, 65 % are less homogeneous based on the drainage and relief 
factor, which implies regions with 25 % that have constant flows in the 
rivers as a result of the high tendency to run-off in the watersheds located 


in the regions of Loreto, San Martin and Cajamarca. 
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Figure 3. Dendrogram obtained by CA for the 91 drainage watersheds 
using Ward Linkage. 


The Drainage and relief are the main cause of greater dissimilarity 
between the studied watersheds. This anomaly is due to the fact that the 
watersheds located on the Atlantic slope and on the Amazon, plains have 
a higher drainage density than the Pacific slope. However, the northern 
area of Peru due to the El Niño Southern Oscillation (ENSO) phenomenon 
is favoured by the intense rains, giving rise to greater floods in the rivers, 
in relation to the southern area where they are scarcer. In previous 


studies by Subyani et al. (2012) and Eltahan et a/. (2021) have obtained 
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in the CA a similarity of the hydrographic watersheds in the geometric 


variables. 


Conclusions 


The prioritization of watersheds is an elemental factor for sustainable 
development and adequate management of water resources. In this 
research, GIS, PCA and WSA techniques were used to characterize the 
morphometric parameters of watersheds with a focus on soil and water 
conservation. The study revealed that the watersheds analyzed in the 
Lima and Amazon regions are the most vulnerable to erosion and present 
potential areas for the application of soil and water conservation. These 
results contribute to decision-making and priority recommendations for 


integrated planning and management of Peru's watersheds. 


The integration of the GIS, PCA and WSA methods was successful. 
The WSA allowed defining the weights for the most relevant parameters. 
Likewise, with the PCA it was possible to determine the largest number of 
significant parameters such as: Lc, Re, Lr, H, K¿and Zmin. 


The homogeneous hydrological characteristics by the dendrogram 
by Ward's method explained that only 35 % of the hydrographic 
watersheds are homogeneous based on geometric parameters and in 
shape, with dissimilarities in sources of water resources. The study 
provides information and useful tools for the development of integrated 


river watersheds management projects at the national level. Also, it 
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provides a series of important data for the implementation of hydraulic 


works for the protection and control of streams in specific regions. 


Abbreviations 


% = Percentage 
AH = Mean stream height 


A = Watersheds area 


C = Constant channel maintenance 

CA = Cluster analysis 

Cov = Covariance of principal component 
CF = Composite factor 

Cv = Coefficient of variation 

PC = Principal component 

PR = Preliminary ranking 

Cr = Roughness coefficient 

D-0,/5 = Matriz diagonal 

Da = Drainage density 

e = Equidistance between contour lines 
ENSO = El Niño Southern Oscillation 
FA = Factorial analysis 


Fs = Form Factor 
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GIS = Geographic Information Systems 
H = Altimetric amplitude 
K¿ = Gravelius compactness coefficient 
km = Kilometer 
Km? = Kilometer squared 
L = Watersheds length 
Lf = Mean length of overland flow 
l; = Contour length 
L, = Stream length 
Max = Maximum 
Min = Minimum 
P = Watersheds perimeter 
PCA = Principal component analysis 
r = Correlation coefficient 


R?2 = Determination coefficient 


Ra = Relief ratio 
R¿ = Circularity ratio 
Re = Elongation ratio 


RP = Rotación ortogonal 
S = Variance 


S2 = Covariance 


Se = Slope of the watersheds 
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Sa = Standard deviation 

S, = Mean stream slope 

Sw = Shape index 

W = Watersheds Amplitude 

WSA = Weighted Sum Approach 

X = Matrix 

XT = Transpose variables 

Z = Vector 

Zmax = Maximum altitude 

Zmin = Minimum altitude 


Z, = Rotated components 
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